On the origins of apparent fragile-to-strong transitions of protein hydration waters 
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2 H NMR is used to study the mechanisms for the reorientation of protein hydration water. In the 
past, crossovers in temperature-dependent correlation times were reported at Ta;i«225K (Xi) and 
T X 2~200K (X2). We show that neither Xi nor X2 are related to a fragile-to-strong transition. Our 
results rule out an existence of Xi. Also, they indicate that water performs thermally activated and 
distorted tetrahedral jumps at T <T X 2, implying that X2 originates in an onset of this motion, which 
may be related to a universal defect diffusion in materials with defined hydrogen-bond networks. 



Studying the temperature dependence of the structural 
relaxation, one can distinguish 'strong' and 'fragile' glass- 
forming liquids, which do and do not show Arrhenius be- 
havior, respectively [lj. For supercooled bulk water, it 
was proposed that a fragile-to-strong transition (FST) 
exists at about 225 K [2j, [3|. However, inevitable crystal- 
lization interferes with direct observation |4| . In confine- 
ment, crystallization can be avoided so that the dynam- 
ics of supercooled water are accessible down to the glass 
transition temperature T g . For confined water, crossovers 
in the T dependence of a correlation time r were observed 
at T xl « 225 K (X x ) and T x2 ~ 200 K (X 2 ) using quasi- 
elastic neutron scattering (QENS) @, Q and dielectric 
spectroscopy (DS) @, HI, respectively. Several workers 
took Xi and X2 as indications for a FST 0, 0. Also. Xi 
was related to a liquid- liquid phase transition |5l. Irj [9I [Io| . 
Challenging these conclusions, others argued that the Ar- 
rhenius processes Pi and P2 observed below Xi and X2, 
respectively, are secondary (/?) relaxations JjJ, , while 
structural (a) relaxation is difficult to observe [8| . 

Confined waters are of enormous importance for bio- 
logical, geological, and technological processes, e.g., an 
interplay of water and protein dynamics enables biolog- 
ical functions [l3j]. Here, we investigate the hydration 
waters of elastin (E) and collagen (C), two proteins of 
the connective tissue. Using 2 H NMR, we study time 
scale and mechanism of water reorientation during T x \ 
and T X 2- The results show that a FST does not exist, 
while tetrahedral jumps become important upon cooling. 

In 2 H NMR, the quadrupolar frequency is probed 14j |: 



± - (3 cos 2 0-1-?? sin 2 9 cos 2(f)) (1) 



Here, the angles (6, <j>) describe the orientation of the elec- 
tric field gradient (EFG) tensor at the nuclear site with 
respect to the external static magnetic field. Since the 
molecular orientation determines the orientation of the 
EFG tensor, rotational jumps render ojq time dependent. 
The anisotropy and asymmetry of the tensorial interac- 
tion are characterized by 5 and 77, respectively. The ± 
signs correspond to two allowed transitions between three 
Zeeman levels of the 2 H nucleus (1=1). 

C, E, and D2O were purchased from Aldrich. Weighed 
amounts of protein and D2O were carefully mixed and 
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FIG. 1: (a) Buildup of the magnetization M z (t) after satura- 
tion for C25 at 230 K (shifted) and C60 at 235 K together with 
fits to two- and three-step relaxations, respectively. The two- 
step fitting function is indicated, (b) T-dependent spectra of 
E43 resulting from the solid-echo sequence, 90j-A-90$J-A-t. 



sealed in the NMR tube to prepare samples with hydra- 
tion levels /i = 0.25— 0.96 (g D2O/I g protein). We refer to 
the samples as 'E' or 'C followed by the value of h in per- 
cent. 2 H NMR spin-lattice relaxation (SLR), line-shape 
(LS), and stimulated-echo (STE) measurements are per- 
formed at a Larmor frequency of ujl = 27rx76.8 MHz. The 
experimental setup is described in Ref. [l5j |. 

In the studied samples, supercoolable bound water and 
freezable free water can coexist. Furthermore, deuterons 
replace exchangeable hydrogens of the protein fl6t . On 
the basis of the amino acid compositions 



17| , one expects 



10-30% of the deuterons to be bound to proteins for the 
used hydration levels. Hence, deuterons of supercooled 
water, crystalline water, and protein can contribute to 
2 H NMR signals. Accordingly, three steps are observed 
in SLR measurements for C60, see Fig. [If a). The third 
step results from deuterons of crystalline water, as can 
be inferred from its absence for low hydration levels, e.g., 
for C25, when free water does not exist amce 
water deuterons outnumber protein deuterons, the higher 
first and lower second steps can be assigned to super- 
cooled water and protein, respectively. The existence of 
distinct steps means that deuteron exchange between the 
three species is slow on the time scale of SLR 0, en- 
abling separate analysis of the dynamical behaviors in 2 H 



2 



NMR. 

Depending on the value of h, we use two- or three-step 
relaxations to fit the buildup of magnetization, M z (t), see 
Fig.QJa). While we will not discuss results for crystalline 
water, Figure^a) compares the T-dependent SLR times 
T\ w and Ti p of deuterons in supercooled waters and pro- 
teins, respectively. T lju shows a very similar minimum for 
all samples, indicating that water dynamics in the hydra- 
tion shells of C and E are highly comparable. Unlike the 
first step, the second step is nonexponential (f3 p « 0.6), 
typical of amorphous solids [13 ] . For C and E, T\ p and, 
thus, the protein dynamics weakly depend on T. 

Using that reorientation of the hydration waters is ba- 
sically isotropic in the vicinity of the minimum, see below, 
Ti w depends on the spectral density according to [14 1 



T^ = (2/15) 5 2 [J(cj l ) + AJ(2cj l )} 



(2) 



Here, J(u>) = J °° i*2(i) cos(ut)dt. The rotational correla- 
tion function (RCF) F 2 (i)oc(P 2 [cos0(O)]P 2 [cos0(f)]) de- 
scribes the time dependence of the Legendre polynomial 
P 2 of the angle 6. For a Debye process, F 2 (t) = exp(— t/r) 
and Jbpp{u) — t/ (1+uj 2 t 2 ). In this case, Eq. ((2]) takes 
the form derived by Bloembergen, Purcell, and Pound 
(BPP) 21]. However, this approach predicts T\ w = 2.4 ms 
at the minimum, at variance with our results in Fig.[^a). 
Thus, a Debye process does not describe the water dy- 
namics, but a distribution of correlation times G(logr) 
exists, as typical of supercooled liquids. To consider such 
distribution, Cole-Davidson (CD) and Cole-Cole (CC) 
spectral densities proved useful in SLR analyses [13 . [22| 



Jcd{u) 



uj 1 sm[PcD arctan(wr C£ ))]/(l+tj 2 Tp £) ) 
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Employing Jcd {Jcc), a width parameter (3cd = 0.22 
{Pec = 0.50) is obtained from T\ w at the minimum. As- 
suming that Pcd {Pec) is independent of T and inserting 
Jcd {Jcc) into Eq. t cd {t cc ) is extracted from T Vw . 
The mean correlation time (tcd) = t~cdPcd [H| of the 
asymmetric CD distribution and the peak position tcc 
of the symmetric CC distribution are shown for E25 and 
E43 in Fig. Etc). At 220-260 K, J CD and J CC yield com- 
parable results, demonstrating insensitivity to the choice 
of the specific spectral density in this range. At lower T, 
non-Arrhenius and Arrhenius (E a — 0.60 eV) behaviors 
are obtained from use of Jcd and Jcc, respectively, so 
that ambiguity about J(uj) hampers analysis. Anyhow, 
neither Jcd nor Jcc yield evidence that Xi exists. 

In the following, we focus on h < 0.25 — 0.43. Then, 
the hydration shells are fully occupied, while freezable 
water hardly exists HI, 2(|. In Fig.QJb), we see that the 
solid-echo spectra of E43 are comprised of two compo- 
nents. Independent of T, the protein component is given 
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FIG. 2: (a) SLR times T\ w and Ti p for deuterons of super- 
cooled waters and proteins, respectively. The dashed horizon- 
tal line is the expectation for the T% w minimum in the case 
of a Debye process, as calculated from Eq. <(2j) for the exper- 
imental values 5 = 2n x 168 kHz and u L = 2n x 76.8 MHz. (b) 
T-dependent integrated intensity S of solid-echo spectra and 
relative height a w of the first SLR step for hydrated E and C 
samples. The lines show the dependence of S on the corre- 
lation time r, as obtained in simulations of isotropic random 
jumps (IRJ) and tetrahedral jumps (TJ), using the experi- 
mental values 5 = 2-k x 168 kHz, 77 = 0.1, and A = 20 /is. 



by a broad spectrum, which is the typical LS in the ab- 
sence of motion 
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For the water component, this LS is 
observed at low T, but the broad spectrum collapses be- 
tween 185 and 215 K resulting in a narrow Lorentzian at 
high T. The narrow spectrum reveals that, at T> 215 K, 
the water molecules show fast (t<1/<J) isotropic jumps 
that average out the anisotropy of the quadrupolar in- 
teraction, while anisotropic motions, e.g., n flips, do not 
result in a Lorentzian and can be excluded. At T < 185 K, 
the broad spectrum indicates that significant water dy- 
namics is absent on the time scale of 1/(5 « 1 /is. Thus, 
Pi, found in QENS works below T x i @, 0], is not the a 
process so that Xi is not a FST. If Pi were the a process, 
a narrow Lorentzian would be observed as 2 H LS down 
to about 140 K due to rpi <C 1/5, see Fig. GUJc), in clear 
contrast to the findings in Fig. [TJb) . The assignment and 
shape of the lines are confirmed when we single out the 
water contribution in partially relaxed (PR) experiments 
(not shown) , in which we do not wait for complete recov- 
ery of M z after saturation, but start acquisition at times 
between first and second SLR steps. 

The total spectral intensity S, determined by integrat- 
ing the solid-echo spectra after correction for the Curie 
factor, exhibits a minimum S « 0.25 at T« 210 K, see Fig. 
|^Jb) , indicating that water dynamics during the dephas- 
ing and rephasing periods A of the solid-echo sequence 
interferes with echo formation [23}. Random- walk sim- 
ulations [23 ] for various motional models show that S 
is a minimum at r = 3 ^s. The signal is not reduced for 
much faster (slower) dynamics, when luq is time-averaged 
(time-invariant). Deviations between measured and cal- 
culated data result because the simulations do not con- 
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FIG. 3: (a) F% c (t m ; t p = 30 (is) of E43 at various T and fits 
to Eq. ©. (b) Fi s {t m \t p = 2, 10,40/Lts) of E30 from PR ex- 
periments at 185 K and fits to Eq. @. (c) Correlation times 
from LS (triangle) and SLR analysis at T>195K (Jcd'- cir- 
cles and squares, Jcc- stars) and from STE analysis below 
195 K [E43: (r cc ) from F 2 cc (t m ; t p = 30 (is), E30: (r„> from 
PR F2 S (t m ',t p = 2(is) ~F2(t m )]. Open symbols mark differing 
SLR results from Jcd and Jcc- Correlation times from mea- 
surements using DS for (solid line) hydrated myoglobin [311 ] 
and (dotted line) E (ft = 0.1)[p|, (+) TSC for E (ft = 0.5) 
[U, (x) MR for C (ft > 0.5) [Tf7and (|) QENS for hydrated 
lysozyme [f|. (d) Evolution-time dependence of (r cc ) (solid 
circles) and (t ss ) (open circles) from PR STE experiments on 
E30 at 185 K and simulation results for (+) distorted (±3°) 
tetrahedral jumps and (*) isotropic 10° jumps. 



sider a distribution G(logr) and protein contributions. 
In Fig. [3jc), we see that LS and SLR analyses yield con- 
sistent correlation times. In SLR experiments, the height 
a w of the step due to supercooled water is also a mini- 
mum at 210 K since a solid echo is used, see Fig.^b). 

Finally, we perform STE experiments [3| to study wa- 
ter dynamics at T < T x %. In STE experiments, we cor- 
relate the instantaneous quadrupolar frequencies luq of 
a deuteron during two short evolution times t p <C r that 
are separated by a longer mixing time t m « r. Using 
appropriate pulse sequences, variation of t m for constant 
t p enables measurement of RCF (£ = sin, cos; x = ss,cc) 



Fg(t m ;t p ) oc (£[uQ(0)t p }Z[ujQ(t m )t p ]) 



(3) 



The brackets (. . .} denote the ensemble avera ge. F* a and 
i^2° result for £ = sin and £ = cos, respectively [l4|. These 
RCF decay when slow (t p w 10 (is < t < T\ « 1 s) molecular 
reorientation alters the value of luq(6, (f>) during t m . 

Figure [3] presents F$(t m ;t p ) of E30 and E43. Water 
dynamics lead to decays at short times ($ w ), while SLR 
results in additional damping of the water (R w ) and pro- 
tein (R p ) contributions at longer times. Taking R w and 



R p from the above SLR analysis, we fit F;f(t m ) to 

A w [(l-B)^ w {t m )+ B}R w {t m ) + {l-A w )Rp{t m ) (4) 

Here, B is introduced to consider that water dynamics 
does not destroy all orientational correlation, see below. 
Using a stretched exponential § w (t m ) = exp[— (tmjr) 13 ], 
we obtain stretching parameters (3 = 0.27 — 0.28 [25| for 
all studied values of t p and T. Thus, hydration water 
exhibits strongly nonexponential RCF below 200 K. Use 
of the r function enables calculation of the mean corre- 
lation time, {t) = (t/(3)T(1/(3). In Fig.^c), we see that 
(t ss ) and (t cc ) from F 2 SS of E30 and F 2 CC of E43, respec- 
tively, follow the same Arrhcnius law (dashed line) with 
activation energy E a = 0.45 eV. We note that F| s of E30 
was obtained in PR experiments to minimize the protein 
contribution, leading to (1 — t4 u) ) = 0.05 — 0.10. 

Comparison of our SLR, LS, and STE data implies a 
crossover in the T dependence of water dynamics in the 
vicinity of T X 2- Here, exact determination of a crossover 
temperature is hampered by a dependence on the choice 
of the spectral density used in SLR analysis. For an 
understanding of the origin of X 2 , knowledge about the 
mechanism for water dynamics below T x i is of particular 
importance. We exploit that the dependence of STE de- 
cays on the length of the evolution time has been shown 
[b€ |2o| to yield valuable insights into the geometry of 
rotational motion since sensitivity to small changes of 
luq and, thus, small angular displacements is higher for 
large t p , see Eq. ([3]). An analysis of (T x (t p )) enables a 
determination of jump angles. While small-angle reori- 
entation, results in a strong decrease of (r x (t c )), (t x ) is 
independent of t p for large-angle reorientation [lj|, see 
Fig.CJd). Extraction of (T x {t p )) from F£{t m ;t p ) of E30 
shows that a substantial dependence on t p is absent, in 
harmony with results for C25 (not shown). Hence, be- 
low T X 2, water exhibits large-angle rather than small- 
angle reorientation typical of the a process [l4[ • In Fig. 
[3Jb), we see that i^ s (t m ; t p ) decays to a small, but finite 
plateau B(t p ), see Eq. ((4]), before the onset of SLR. The 
plateau height depends on the geometry of the motion, 
e.g., B(t p — >oo) = 1/n for a n-site jump [26|. Here, B de- 
creases from 0.16 to 0.06 when t p is extended from 2 to 
40 (is, indicating that, though water reorientation is not 
isotropic, angular restrictions are not severe j26|. While 
exact two- (B=l/2) or four-site (.6 = 1/4) jumps can be 
excluded, distorted tetrahedral jumps, which may be ex- 
pected for a disordered hydrogen-bond network, are con- 
sistent with the data In 2 H STE work on ice I h [27j], 
a similar dependence B(t p ) was shown to indicate that 
translational diffusion of water molecules via interstitial 
defects involves distorted (±3°) large-angle reorientation 
between 7 O-D bond orientations in the crystal. 

In summary, we exploited the capabilities of 2 H NMR 
to ascertain correlation times and mechanisms for the ro- 
tational motion of supercooled protein hydration waters. 
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In the literature @, @, 0, @, fl EH EL El , it was a contro- 
versial issue to take crossovers Xi and X2 of temperature- 
dependent t as evidence for a FST. Prerequisite for such 
interpretation of Xi or X 2 is that the processes Pi or P 2 
below these crossovers are the a process. The present 
results, e.g., the LS at 145—185 K, rule out that Pi is the 
a process and, hence, that Xi at T x \ m 225 K is a FST. 
While a narrow spectrum would result if Pi were the a 
process, a broad spectrum is observed. Our fin ding s con- 
flict with QENS (E i] and X H NMR diffusion Q stud- 
ies, attributing Xi to a FST. Further evidence against Xi 
being a FST comes from extrapolation of rp% (T) , which 
yields too small values T g < 100 K. Also, a weak wave- 
vector dependence t(Q) below Xi |6j implies local rather 
than diffusive motion. Our results are consistent with 
the conjecture that Pi is a (3 process [1, 11, 12|, if the 
underlying motion is too restricted to be probed in 2 H 
NMR. 

Exploiting the sensitivity of 2 H NMR to the mecha- 
nisms for water reorientation, we investigated the origin 
of X 2 at T x2 ~20QK, reported in DS work dED]. Above 
T X 2, SLR and LS analyses showed that water exhibits 
isotropic reorientation described by a broad distribution 
G(logr). Below T X 2, STE experiments indicate that wa- 
ter performs large-angle jumps, most probably distorted 
tetrahedral jumps, which follow an Arrhenius law with 
E a = 0.45 eV. For P 2 , DS work on water in various 
confinements [ll| reported comparable rp 2 (T), see Fig. 
[3jc). Also, mechanical- relaxation (MR) and thermally- 
stimulated currents (TSC) studies found a process with 
similar values of r [lil l2(ij|. Therefore, all these methods 
probe P 2 . However, P 2 neither destroys all orientational 
correlation, see Fig. [31 nor affects rigidity Thus, 
P 2 is not the a process and X 2 is not a FST. More- 
over, the large-angle jump mechanism is strong evidence 
against P 2 being a Johari-Goldstein j3 process [29( , which 
results from small- amplitude reorientation [3fjl | . Compar- 
ison with previous 2 H NMR work [27j showed that diffu- 
sion of water molecules via interstitial defects is consis- 
tent with the present results. Since E a «0.45eV is not 



only found for water in various confinements 11| , but also 
for crystalline and glassy bulk water Q , we suggest that 
water shows a characteristic tetrahedral jump motion, 
which is controlled by breaking of hydrogen bonds and, 
possibly, related to interstitial defect diffusion whenever 
a defined hydrogen-bond network is established, although 
we cannot exclude that the tetrahedral jump motion is 
governed by the protein surfaces in our case. The onset 
of this tetrahedral jump motion leads to more or less pro- 
nounced crossovers X 2 . Extrapolating rp 2 (T), a value of 
100 s is reached in the vicinity of 136 K, the first widely 
accepted [3| , but later questioned Q value of T g for bulk 
water. Thus, one might speculate that the reported small 
calorimetric effects are not related to a glass transition, 
but to freezing of interstitial defects. 

2 H NMR, which probes single-particle RCF, cannot 



resolve the issue whether the a process exists below X 2 , 
because P 2 destroys basically all orientational correlation 
before an onset of structural relaxation. Detection of 
the a process rather requires techniques sensitive to the 
reorganisation of the whole hydrogen-bond network. In 
this respect, it is interesting that MR and TSC studies 



20j found a slower process P3, see Fig. [3jc), which 
affects rigidity and, hence, may be the a process. Finally, 
it is an open question whether the onset of tetrahedral 
jump motion is related to a liquid-liquid phase transition, 
which was proposed to lead to a low-density liquid with 
a more ordered tetrahedral network upon cooling 0, Q • 

Funding of the DFG through grants VO 905/3-1 and 
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